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We address the problem of finding self-consistent parquet equations for two-body corre- 
lated Greens functions that arise out of a cumulant expansion. The general theory as de- 
veloped for non-relativistic electron systems reproduces the conventional two-body parquet 
diagrammatic technique for the vertex functions. As an application, the crossing-invariant 
linearized basic parquet equations are treated for the correlated two-particle Greens function. 
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§1. Introduction 

Parquet equations are used to self-consistently sum up ring and ladder diagrams 
that are needed for the microscopic description of numerous many-body systems 
including the Coulomb gas, liquid helium, and nuclear matter. These equations 
were originally derived by Landau et al. 1 ) to study the nonperturbative domain in 
quantum electrodynamics and in the theory of two-particle scattering to construct an 
extension of the Bethe-Salpeter equation 2 ) that accounts for the crossing symmetry, 
which is a consequence of the Pauli principle. This symmetry is easily described, but 
extremely difficult to implement in practical applications. For a two-particle Greens 
function (GFs) that refers to Fermions, this fundamental symmetry is expressed by 

G(12, 1'2') = -G(21, 1'2') = -G(12, 2'l') = (7(21, 2'l'), (M) 

where the numbers collect space, imaginary time, and spin indices (e.g., 1 stands 
for the space r%, time t±, and spin s\ coordinates). It is important that this simple 
symmetry couples different scattering channels to each other by means of partial irre- 
ducible vertex functions. This coupling is accounted for by a closed set of nonlinear 
self-consistent equations for the two-particle irreducible vertices, which are called 
parquet equations. The theory is applicable not only for small but also for large 
values of the coupling constants and allows a simultaneous treatment of long-range 
as well as short-range correlations. 

Applications of the parquet approach start from an initial set of irreducible di- 
agrams, which cannot be broken into two parts by cutting a pair of electron lines to 
generate all related reducible diagrams of the respective n-body GF. The full par- 
quet equations, which are based on all irreducible diagrams, are exact, in accordance 
with the crossing symmetry, and conserving with respect to the particle number, 
the energy, the momentum, and the angular momentum. Unfortunately, this in- 
ternal consistency of the full parquet theory breaks down in any approximation. 7 ) 
Since the rigorous solution of the full equations is not feasible for even the simplest 
models, it is necessary to treat simplified, manageable variants of the originally ex- 
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act full theory. Several diagrammatic approximations beyond perturbation theory 
have been proposed in the literature and successfully applied to the treatment of 
many-body problems. 8 )~ 13 ) A further systematic extension of all these techniques is 
possible by a calculational scheme, which is based on a rearrangement of the whole 
diagrammatic collection. Such a transformation of the entire hierarchy of equations 
for the many-body GFs is carried out by using an alternative generating functional. 
Of particular interest in statistical physics is the cumulant expansion, which pro- 
vides a natural framework for the definition of an n-particle approximation as being 
the set of equations for correlation functions that terminates, when cumulants of 
the order n + 1 are discarded. Applied to the many-particle theory, the cumulant 
expansion generates equations of motion for correlated thermodynamic GFs that is 
nothing but a rearrangement of the Kadanoff-Baym hierarchy. 14 ) It is the advantage 
of this cumulant approach that it provides a firm basis for a systematic self-consistent 
n-particle approach by cutting the chain of equations at the (n + l)-particle level. 
This rearranged diagrammatic expansion has to satisfy the crossing symmetry, which 
is a fundamental property of any many-particle system. Unfortunately, the result- 
ing set of equations violates the crossing symmetry. This observation motivates 
our work to propose a system of crossing-invariant parquet equations for correlated 
thermodynamic GFs that are particularly useful for the construction of n-particle 
approximations . 



§2. Correlated Green functions 



Let us develop the many-body theory for an electron gas in a solid that is 
described by the Hamiltonian 

H = Y.j dr 4(r,t) -^-A r + V{r) 

S L J 

drdr'^tir^^ir^tMr-r'^ir'^Mr^), (2-1) 

s,s' 

with ipt(r,t) [ifj s (r,t)] denoting the creation [annihilation] operator for a Fermion 
with spin s that occupies the lattice site r at time t. The carriers with an effective 
mass m move in the single-particle potential V{r) and interact with each other via 
the Coulomb potential v(r — r'). The equal-time anticommutator relations of the 
field operators have the form 



+ 



r,t),^ s ,(r',t) 
$ s ,s>5(r - r'), 



(2-2a) 
(2-2b) 



while the time dependence of the system is governed by the von Neumann equation 



- ih—il> a (r,t) = [H,ip s (r,t)]_ , 



(2-3) 
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which is used to derive equations of motion for the GFs. A straightforward approach 
is based on the generating functional 

oo 

G[A, rj]=l + Y^ 1-^2 / dl • ' ' dnd1 ' ' ' ' dn ' A ( n ) • • • A (!) 
n =i ^ n '> •* 

xG(l . . . n, l' . . . n')rj(l') . . . i](n'), (2-4) 

in which both r](j) and X(j) denote anticommutating functions. The infinite chain 
of coupled equations for the n-particle GFs G(l ... n, 1' ... n') is generated from a 
functional differential equation for G[X, rj] 

^-•w^}^^ (2 ' 5) 

with V(l — 2) = v(ri — r2)6(t\ — £2)- An alternative hierarchy of equations is 
derived for correlated n-body GFs G c (l . . . n, V . . . n') by introducing its generating 
functional via the equation 

G[X,r,]=exp(G c [X,r,]). (2-6) 

Exploiting Eq. (|2-5p , the exact results for the one-particle and two-particle correlated 
GFs are easily derived. 14 ) The first equations of this infinite chain have the form 

G(l, 1') = G (l, 1') -ihj dId2G (l,lMT - 2)G C (2T, 2 + l'), (2-7) 



G c (12, 1'2') = ih J dld2i/(l-2)j-G c (212,2 + l'2') 
+ [Go(l,T)G(2,2)G(T,l')G(2,2') - G (1,T)G(2,2)G(T,2')G(2, l' 



+ 



(2-8a) 

')] (2-8b) 
(2-8c) 

(2-8d) 
(2-8e) 

where the appearance of the three-particle GF in Eq. (|2-8aj) should be noted. It is a 
special feature of this set of equations for the cumulants that the Hartree-Fock GF 

Go appears in each equation of the hierarchy ''• ; " " " ' " f ' '•'•■■• u ' '•■ 

determined by a integro-differential equation 



+ [G (l,l)G(2,2)G c (12,l'2, j 

G (1,T)G(2, l')G c (T2, 2 + 2') + G (1,T)G(2, 2')G C (12, l'2+) 

G (1,T)G(T, l')G c (22, 2 + 2') + G (1,T)G(T, 2')G C (22, l'2+) 



^A ri - V(n) + ihf &1V(1 - T)G(T, 1 + ) + ih^ G (l, 1') 
-ih J dTF(l-T)G(l,T + )G (T, 1') = 6(1-1'), 



(2-9) 
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which includes the full one-particle GF G in a self-consistent manner. Taking into 
account the relationship 

G(12, 1'2') = G c (12, 1'2') + G(l, l')G(2, 2') - G(l, 2')G(2, l') (2-10) 

between the standard and correlated two-particle GF, we conclude that the effective 
one-particle approach as given by the Hartree-Fock approximation results, when the 
correlated two-particle propagator G c (12, 1'2') is neglected. Accordingly, an opti- 
mized pair approximation is established, when the three-particle GF in Eq. ()2-8al) 
is removed. The resulting closed equation for the correlated two-particle GF dif- 
fers from the related standard result by the appearance of the GF Go- Indeed, by 
replacing Go by G in Eqs. (|2-8a|) to (|2-8ej) . we recover the linearized basic parquet 
equation, in which only the lowest-order irreducible vertex function is taken into 
account. Besides this formal correspondence, the crossing symmetry is basically 
violated in Eq. (2.8). This serious defect of the cumulant expansion has to be re- 
moved before this powerful approach can be exploited. What is needed is a crossing 
invariant theory for correlated many-particle GFs. 

§3. Parquet equations 

The full imaginary-time parquet equations self-consistently generate all reducible 
Feynman diagrams for the two-body GF in terms of the bare interaction, the fully 
dressed one-particle GF, and an initial set of irreducible diagrams by maintaining the 
complete crossing symmetry. The crossing operation itself transfers reducible vertex 
diagrams of a given scattering channel into irreducible vertices of complementary 
channels so that a coupled, non-linear set of self-consistent integral equations arises 
for the irreducible two-particle vertices. To apply this calculational procedure to 
the hierarchy of correlated GFs, we chose a mixed representation, in which not only 
vertex functions appear, but also partial two-particle GFs. Accordingly, the expres- 
sion for the total correlated four-point GF G c decomposes into a sum of different 
scattering contributions 

G c = go + g ee + 9d + 9c- (3-1) 

The collection of all irreducible vertex diagrams I] rr leads to the component go given 
by 

50 (12,l'2 / ) = ~y dT...d4yl (12,r2)r irr (T2,34)G(3,l')G(4,2 / ), (3-2) 

where Aq denotes the symmetrized product of one-particle GFs 

4,(12, 1'2') = G (l, l')G(2, 2') + G(l, l')G (2, 2'). (3-3) 

Particle-particle ladder diagrams are summed in g ee , which satisfies the self-consistent 
integral equation 

ffee (12,l'2') = i I dT...d4ylo(12,T2)r ee (T2,34)G c (34,l / 2'). (3-4) 
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The construction of the particle- hole hierarchy is more involved. To account for the 
crossing symmetry, we make use of the decomposition 

^ = 5 (3-5) 
where the components satisfy the respective integral equations 

5 i 1) (12,l'2 / ) =- J dl...d4G (l,T)G(3,l , )r (i (T2,34)G c (42,22'), (3-6) 
5 «(12,l'2 / ) = ydT...d4Go(l,T)G(4,2 / )r c (T2,34)G c (32,l / 2). (3-7) 

(2) 

The quantities g d c are obtained from the right-hand side of these equations by 
interchanging the numbers 1 with 2 as well as 1' with 2'. The vertex functions, 
which enter Eqs. (13-4p . (|3-6p . and (!3-7p . are given by 



r ee = r irr + 7 d + 7c, r d = r irr + 7 C + 7 ee , r c = r irr + 7 d + j ee , (3-8) 

where the partial vertices j a (with a: = c, cZ, or ee) generate the GFs g a via the 
relationship 

5a (12,l'2 , ) = y"dT...d4G(l,T)G(2,2) 7a (T2,34)G(3,l')G(4,2'). (3-9) 

The set of Eqs. {HD} to (EZHJ) together with Eqs. fl22j) and ([231) establish the self- 
consistent parquet approach for the calculation of the correlated one- and two- 
particle GFs. An iteration of these equations can start from an effective expression 
for the Hartree-Fock GF Go and an initial set of irreducible vertex diag rams -Z~irr that 
is used to calculate the first approximation for G c from Eqs. (|3-ip and f|3-2[) . To pro- 
ceed further, also the partial vertex functions 7^, 7 C , and 7 ee have to be determined 
from the integral Eq. (|3-9|) . The solution of these equations is facilitated by a Fourier 
transformation with respect to the spatial coordinates and by a transformation to 
Matsubara frequencies z v according to 

G c (12,l'2')= {J-j^j G c (xiX2iX V x 2 '\z v z v iU) n ) (3-10) 



v,V ,n 



x exp \-%z v t\ - i{uj n - z u )t 2 + iz v dy + i{uj n - z v >)t2>] , 

where j3 = 1/(/cbT) denotes the temperature parameter and z v = mv/(h(3) + [i/h, 
uj n = iirn/(hP)+2fj,/h, with \x being the chemical potential. The result is an algebraic 
relationship 

9s lS 2,s3s 4 ( kk 'q\z v z v iuj n ) = G(k,z v )G{q -k,u n - z v ) 
xis^i^sAkk'qlzvZv'UJ^Gik' , z v i)G(q - k',uj n - z v >), (3-11) 

which is easily inverted. Furthermore, the 7 C diagrams are obtained from the ones 
for 7^ by exploiting the antisymmetry of the vertex function under the exchange of 
the electron (or hole) variables. This specific crossing symmetry is expressed by 

7 d (12,l'2 , ) = - 7c (12,2 , l'). (3-12) 
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The set of parquet Eqs. (j3-ip to (|3-9p for the correlated two-body GF is the main 
result of this paper. The solution of these equations together with Eqs. (|2-7|) and (|2 9j) 
is a formidable task that is generally not possible without further approximations. 
Therefore, we discuss two limiting cases, which have received particular attention. 

First, let us derive the conventional parquet equations for GFs that have been 
studied in a number of review articles. 3 ^ 6 ^ As already mentioned, the main pecu- 
liarity in the theory of correlated GFs is the appearance of Go- As the linearized 
basic parquet equations are obtained from Eq. (2.8) by the replacement Go — ► G, 
one expects a similar unfolding for the set of Eqs. (|3-ip to (|3-9p . Indeed, replacing 
Go by G in these parquet equations, we recover closed equations for the total vertex 
function r, which is defined by Eq. (|3 9j) . with g a and j a being replaced by G c and 
r, respectively. Obvious manipulations lead to the final equations 

r = r irr + 7 ee + 7 d + 7 c , (3-13) 
7d (12, 1'2') = - J dT...d4T d (lT, l'2)G(2,3)G(4,T)r(32,42'), (3-14) 
7 C (12,1'2') = ydT...d4r c (lT,22')G(2,3)G(4,T)r(32,l'4), (3-15) 

7 ee (12,l , 2 / ) = ±J dT...d4r ee (12,T2)G(T,3)G(2,4)r(34,l'2 / ), (3-16) 

which agree with the conventional parquet equations for the vertex function. 8 ) Again, 
the vertices r c , i - ^, and F ee of different scattering channels appear, which are defined 
by Eq. (|3-8p . The exact correspondence between the parquet equations for the two 
different sets of many-body GFs (namely G and G c ) provides strong evidence for 
our crossing invariant result in Eqs. (|3T|) to (|3-9|) . 

For practical applications, a linearization of the full parquet equations is some- 
times desirable. A further simplification consists in a restriction to the lowest-order 
irreducible vertex function, which is given by 

r (12, 1'2') = ihV(l - 2) [5(1 - l')5{2 - 2') - 5(1 - 2')5(2 - l')] . (3-17) 



Adopting both approximations, the parquet Eqs. (|3-ip to (|3-9p lead to an equation 
for the correlated two-particle GF, which is diagrammatically represented in Fig. 1. 
The collection of all diagrams in the rows I to IV satisfies the crossing symmetry 
and has, therefore, the potential to describe distinct physical entities. Particle- 
particle and particle-hole ladder diagrams are generated by the channels II and (III, 
IV), respectively. A mixing of spin states appears in the electron-hole channel IV, 
which encompasses components that originate from the vertex functions 7^ and 7 C . 
The equation for G c as represented by the diagrams in Fig. 1 is nothing but a 
crossing invariant extension of Eq. (2.8), which was straightforwardly derived from 
the generating functional G c [X,rj\. 

To further discuss the physical significance of our result, let us treat two examples 
in more detail. The set of diagrams shown in Fig. 1 by row I fulfill the crossing 
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Fig. 1. Feynman diagrams for the correlated two-particle GF G c (12, 1'2') arranged in four rows 
I to IV, which collect subgroups that satisfy the crossing symmetry. The bare potential V is 
represented by a dashed line. One-particle propagators G and Go are shown by oriented solid 
lines. 

symmetry. All contributions are composed from one-particle GFs that are combined 
in the GW approximation for the self-energy. In the Fourier space, we have 

E(k,z v ) =G 1 (k,z u )-G- 1 (k,z u ) = -J2w{k,q\cj n )G(q - k,u n - z„), (3-18) 

<j,n 



w(k, q\u n ) =-p H fc ~ fcl ) 2 + v ( k + kl ~ q ^ 2 ~ v ( k ~ k i) v ( k + k i- q)] 



hi, v\ 



xGo(fci, z Ul )G(q - ki 



(3-19) 



This result can be viewed as a natural two-particle extension of the Hartree-Fock 
approximation (G = Go), which trivially satisfies the crossing symmetry. The full 
self-consistent GW approximation in Eqs. (|3- 18|) and (|3- 19|) has the special feature 
that the kernel W depends on the squared coupling constant as well as on the 
Hartree-Fock GF G . 

More interesting is the crossing symmetric Cooper channel II, which collects all 
relevant diagrams of the BCS theory of superconductivity. As the V dependent parts 
of 7 ee are never large, 8 ) we restrict ourselves to the homogeneous integral equation, 
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which is expressed by 

G c (12, 1'2') =iKj dId2F(T - 2) (3-20) 
xi [G (1,T)G(2, 2) + G(1,T)G (2, 2)] G c (12, l'2'). 



An analytic solution of this equation together with Eq. (|2-7j) is possible by exploiting 
the symmetry relation 

G c (xiX2,xi'X 2 '|^^'u; n ) = -G*(xiix 2 t,xiX2\ztiZ*u*), (3-21) 

(with Xj = rj, Sj) and by introducing the gap function A according to 

A(xxx 2 ,x v x 2 i\u n ) =v(r 1 -r 2 )v{r v - r 2l )^ s ^G c {x\x 2 ,x v x 2 >\z v z l ,iu ri ). (3-22) 

Applying a spatial Fourier transformation, the coupled Eqs. (I2-7|) and (13-2(1 are 
expressed by 

G(k, z v ) = G (k, z u ) - G (k, z u )^- e £ ^A(kk'q\oj n ) (3-23) 

q,n 

[G (k,z u )G(q - k,uj n - z v ) + G(k,z u )G (q -k,u n - z v )\ , 

A{kk'q\uj n ) = ^2v(k- kx)A(k-yk'q\u n ) (3-24) 
fci 

X {~2^ [^0(^1,^)^(9 - k\,u n - z v ) + G(ki,z v )G (q - kx,u n - z u ) 

with e — » +0. These equations for the one-particle and gap functions provide the 
basis for the study of BCS suerconductivity. To simplify the analysis further and to 
derive approximate analytical results, it is assumed that the Cooperons exist only in 
a very narrow zero-frequency mode so that the gap function is given by 

A(kk'q\u n ) = 2 XI A Sl s 2 ,s lS2 (kk'q\uj n ) = -ih(35 nfl 5 q pA(k,k'). (3-25) 

81,82 

Adopting an effective expression for the Hartree-Fock GF of the form 

with e(k) denoting the kinetic energy of carriers and restricting to a short-range 
attraction with the coupling constant g, the gap equation reduces to 
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which is exactly the result of the well-known BCS theory. 

Recapitulating this derivation, it is tempting to speculate that in the crossing 
symmetric particle-hole channels III and IV there also exist specific, physically mean- 
ingful two-particle excitations such as, e.g., the paramagnon. 8 ) Both a numerical and 
analytical study of all these contributions to the correlated two-particle GF, which 
are given by 

G c (12, 1'2') ~ihj dTd2U(I - 2) (3-28) 
xi |g (1,I)G(2, l')G c (T2, 2+2') + G (2,T)G(2, l')G c (lT, 2 + 2') 
+G (1,T)G(2, 2')G C (12, l'2 + ) + G (2,T)G(2, 2')G c (lT, l'2 + ) 



G c (12, V2') ~ -ih J dId2U(T - 2) (3-29) 
x1|g (1,T)G(T, l')G c (22, 2 + 2') + G (1,T)G(T, 2')G C (22, l'2 + ) 
+G (2,T)G(T, l')G c (12, 2+2') + G (2,T)G(T, 2')G C (12, l'2 + ) 



for the particle-hole channels, would contribute to a systematic analysis of the pair 
approximation in the many-body theory. In particular, such an approach would be 
applicable to critical regions of phase transitions in the vicinity of singularities in 
correlation functions and would allow the identification of symmetry breaking order 
parameters. 



§4. Summary 



The conventional parquet equations couple two-particle scattering channels, which 
are described by partial vertex functions of thermodynamic GFs to account for the 
crossing symmetry. This many-body theory is completely exact, when all irreducible 
vertices are included in the irreducible vertex function. Unfortunately, the full par- 
quet equations cannot be solved without truncations. Any approximate treatment 
of the parquet equations, however, depends sensitively on the grouping of correlation 
functions in the hierarchy of equations. An attractive alternative to the conventional 
equation of motions for thermodynamic GFs represents the cumulant expansion, 
which generates a chain of equations for correlated GFs. It is an advantage of the 
G c hierarchy that it creates step by step self-consistent n-particle approximations. 
For instance, the pair approximation arises in a natural way, when the correlated 
three-particle GF is neglected. Unfortunately, this procedure has a serious defect: 
the resulting equations for the correlated GFs do not exhibit crossing symmetry. 
Therefore, an extension of basic equations for correlated GFs is needed that is in 
line with this fundamental symmetry. It was the aim of this paper to present a 
set of parquet equations that represents a crossing-invariant cumulant expansion, 
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from which the conventional parquet equations are straightforwardly recovered. Of 
particular interest is the linearized version of this parquet approach for correlated 
GFs. 
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